T-Lyrl-17236: A Long-Period Low-Mass Eclipsing Binary 



Jonathan Devor^'^, David Charbonneau^'^, Guillermo Torres^, 
Cullen H. Blake^, Russel J. White^, Markus Rabus^, Francis T. O'Donovan^, 
Georgi Mandushev^, Caspar Bakos^, Gabor Furesz^, and Andrew Szentgyorgyi^ 

66" 

O ■ ABSTRACT 
O 
(N 

We describe the discovery of a 0.68+0.52 Mq echpsing binary (EE) with an 
8.4-day orbital period, found through a systematic search of ten fields of the 
Trans-atlantic Exoplanet Survey (TrES). Such long-period low-mass EEs consti- 
tute critical test cases for resolving the long standing discrepancy between the 
^ . theoretical and observational mass-radius relations at the bottom of the main 

sequence. It has been suggested that this discrepancy may be related to strong 
stellar magnetic fields, which are not properly accounted for in current theoreti- 
. cal models. All previously well-characterized low-mass main sequence EEs have 

periods of a few days or less, and their components are therefore expected to 
be rotating rapidly as a result of tidal synchronization, thus generating strong 
^ , magnetic fields. In contrast, the binary system described here has a period that 

Tj- , is over three times longer than previously characterized low-mass main sequence 

EEs, and its components rotate relatively slowly. It is therefore expected to 



(N 
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rn _ 

[-^ , have a weaker magnetic field and to better match the assumptions of theoretical 

O \ stellar models. Our follow-up observations of this EE yield preliminary stellar 

O ' properties that suggest it is indeed consistent with current models. If further 
observations confirm a low level of activity in this system, these determinations 

^ . would provide support for the hypothesis that the mass-radius discrepancy is at 

^ ! least partly due to magnetic activity. 
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Introduction 



Despite a great deal of work that has be en done to understand th e structure of low 



mass (< 0.8 Mq) main sequence stars (e.g., IChabrier fc Ba.raffel 120001). models continue 



to un derestimate their radii by as much as 15% (ILacy 



1977a 



Torres fc Ribas 2002; Ribas 



20061 ) . This is a significant discrepancy, consider ing that for solar-type stars the agreement 
with the observations is typically within 1-2% (lAndersenl Il99ll . Il998l ). In recent years an 
intriguing hypothesis has been put forward, suggesting that strong magnetic fields may have 



bloated these stars, eithe r through chromosph eric activity (e.g., iRibad l2006l : [Torres et al. 



20061 : iLopez-MoralesI 



2007 



disruption (ITorres et al. 



Chabrier et al.ll2007l ) or through magnetically induced convective 



2006h . Such strong magnetic fields are expected to be formed by the 



dynamo mechanism of rapidly rotating starsE To test this hypothesis, one needs to measure 
both the masses and radii of low-mass stars, which thus far can be done most accurately with 
eclipsing binary (EB) systems. However, all well characterized low-mass main sequence EEs 
have orbital periods shorter than three days (see T able [T| and are th erefore expected to have 
synchronization timescales shorter than ~100 Myr (jZahnlll977l . Il994l . see Figure[T]and further 
description in §[6]). As a result of these short periods and short synchronization timescales, 
the rotations of these binary components are expected to have accelerated to the point that 
they now match the rapid angular velocity of their orbits. With such rapid rotations, these 
binary components could have a wide range of dynamo-induced magnetic field strengths. 
To better constrain current stellar models, we set out to find systems with slowly rotating 
components. Such systems would presumably have comparably weak magnetic fields, thus 
being more consistent with the model assumptions. Furthermore, by comparing the mass- 
radius relations of binary components with well determined levels of magnetic activity, one 
could test various magnetic disruption models. 

We note here that in addition to EB analysis, long-baseline optical i nterferometry has 



also been used recent l y to measure the radii of nearby low-mass stars (ILane et al.l 12001 



Segransan et al.l l2003l : iBerger et al.l l2006l ) . While these stars are single and are therefore 



^Dynamo theory predicts that this mechanism operates only in partiaUy convective stars. However, 
the strong magnetic activity observed in fully conve ctive low-mass stars indi cates that they also possess a 
mechanism for generating strong magnetic fields (see iBrowning fc Basrill2007l . and references therein). 
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expected to rotate slowly, their masses can only be estimated through empirical mass- 
luminosity relations or other indirect methods. Those determinations are thus less fun- 
damental, in a sense, and arguably of lesser value for accurately constraining stellar models 
and testing the magnetic disruption hypothesis. 



Initial Photometric Observations 



T-Lyr 1-17236 was first identified as a likely low-mass EB candidate in the iDevor et al. 



(120081 ) catalog, following a systematic anal ysis of the light cur ves (LCs) within ten fields of 
the Trans-atlantic Exoplanet Survey (TrES: lAlonso et al.ll2004l ). TrES employs a network of 
three automated telescopes to survey 6° x 6° fields of view. To avoid potential systematic 
noise, we performed our initial search using data from only one telescope , Sleuth, located at 
the Palomar Observatory in Southern California (10 'Donovan et al.ll2004j ). and we combined 
additional data at subsequent follow-up stages. Sleuth has a 10-cm physical aperture and 
a photometric aperture radius of 30". The number of LCs in each field ranges from 10,405 
to 26,495, for a total of 185,445 LCs. The LCs consist of ~2000 Sloan r-band photometric 
measurements binned to a 9-minute cadence. The calibration of the TrES images, identifica- 
tion of stars therein, and the extraction and decorrelation of the LCs are des cribed elsewhere 
dPunham et aPbooi iMandushev et allbood : lo 'Donovan et ahlbood . I2OO7I 1 



An automated pipeline was used to identify and characterize th e EEs among the T rES 
LCs. This pipeline has been described in detail in a previous paper (jPevor et al.ll2008l ). At 
the heart of this anal ysis lie two computational tools: the Detached Eclipsing Binary Light 



curve 



fit teijl (^DEBiL; Devor 2005), and the Method for Eclipsing Component Identificatioiifl 



(MECI; bevor fc Charbonneaul l2006al jbll DEBiL fits each LC to a geometric model of a 
detached EB that consists of two luminous, limb-darkened spheres that describe a Newtonian 
two-body orbit. MECI th en incorporated some of the DEBiL results, and together with 
2MASS color information (ISkrutskie et al.ll2006l ). refit each LC to a physi cal model that 



is constrained b y the solar metallicity Yonsei-Yale theoretical isochrones (lYi et al.l 12001 



Kim et al.ll2002l ). Thus, using only photometric data, the DEBiL/MECI pipeline provided 
initial estimates of the absolute physical properties of each EB. These estimates were then 
used to locate promising candidates for follow-up. 



^The DEBiL source code, utilities, and example files are available online at: 
http://www.cfa.harvard.edu/~jdevor/DEBiL.litml 

^The MECI source code and running examples are available online at: 
[http://www.cfa.harvard.edu/~jdevor/MECLhtml^ 
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Using this pipeline a total of 773 EBs were identified within the TrES dataset. Of these, 
427 EBs were both detached and had small out-of-eclipse distortions, thereby enabling the 
DEBiL/MECI pipeline to estimate their component masses. These results, together with 
many other properties, are listed for each EB in an online catalog f Pevor et al. 2008 ). Of 
these characterized EBs, we then identified a handful of promising long-period low-mass can- 
didates and chose one, T-Lyrl-17236 (aaooo = 19'^07™16^621, (52ooo = +46°39'53".21, P = 
8.429441 ± 0.000033 days; see Table [2] for additional information), for further follow-up 
and an alysis. As with all o f our low-mass candidates, we repeated the MECI analysis us- 
ing the iBaraffe et al.l (119981 ) solar-metallicity isochrones (with a mixing length parameter of 
ttML = 1-0), which are more accurate than the Yonsei-Yale isochrones in this regime. The 
resulting MECI mass-mass likelihood contour plot of T-Lyrl-17236 is shown in Figure [2j 
Since the MECI analysis incorporates data from theoretical stellar models, we cannot use it 
to constrain stellar models. Rather, once we identified the candidate, we followed it up pho- 
tometrically and spectroscopically, and used only these follow-up data to derive the binary's 
absolute properties. 



3. Follow-up Photometric Observations 



In order to characterize T- Lyrl-17236 we comb ined photometric data from four tele- 
scopes: (1) Sleuth and (2) PSST (jPunham et al.ll2004l) of the TrES n etwork, (3) the Instituto 
de Astroffsica de Canarias telescope fIAC80: Galan fc Coboj 1987 ). and (4) the Hungarian 
Automated Telescope Network (HATNet; iBakos et al.ll2004l ). With the exception of the 
IAC80, we obtained our photometric data from archived survey datasets that were intended 
for locating exoplanets. 

As part of the TrES network (see §[2]), Sleuth and PSST are operated similarly. However, 
PSST, which is located at the Lowell Observatory in Arizona, observes in the Johnson i?-band 
whereas Sleuth observes in the Sloan r-band (see Figures [3] and H]). Furthermore, PSST has a 
20" photometric aperture radius compared to Sleuth's 30" radius, which provides PSST with 
a higher resolving power than Sleuth. However, the smaller aperture of PSST also causes it 
to have noisier photometry, with an RMS of 0.031 mag for T-Lyrl-17236, compared to the 
Sleuth photometry that has an RMS of 0.028 mag. Though these differences are small, they 
would have affected our analysis. We therefore chose not to use the PSST data for fitting 
the photometric model, though we did use them to improve the determination of the orbital 
period and the epoch of eclipse (see §[5]). 



"^http: / /www. cfa. harvard. edu/'--^jdevor/Catalog. html 
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In an effort to better constrain the eclipses of T-Lyr 1-17236, we obtained data from tlie 
IAC80, an 82-cm aperture telescope with a 14' x 14' field of view, located at the Observatorio 
del Teide in the Canary Islands. We produced an /-band LC at a 1.3-minute cadence 
using the 1024 x 1024-pixel Tromso CCD Photometer (TCP), resulting in 0.008 mag RMS 
photometry for T-Lyrl-17236. Unfortunately, we were only able to observe a primary eclipse 
with the IAC80. We therefore incorporated archival HATNet observations so as to provide 
coverage of the secondary eclipse in a similar bandpass (see Figures [3] and [5]) . 

HATNet is a network of six 11-cm aperture, fully- automated telescopes (HATs) located 
at the F. L. Whipple Observatory in Arizona and at the Submillimeter Array site atop Mauna 
Kea, Hawaii. The HATs have an 8° x 8° field of view, a response that peaks in the J-band, 
and operate at a 5.5-minute cadence. To reduce the photometric noise, the HAT point spread 



function (PSF) is broadened to a ~15" aperture radius through microstepping (IBakos et al. 



2OO2I ). Even so, the HATNet photometric RMS for T-Lyrl-17236 was comparably large, at 



0.084 mag. Nevertheless, to provide more complete coverage of the primary and secondary 
eclipses in the /-band, we combined the IAC80 observations with data from HAT-7 (Whipple 
Observatory) and from HAT-8 (Mauna Kea). Due to the very different characteristics of 
these two systems, however, we chose not to adopt any of the model parameters derived 
from these data, and only used these results as an independent confirmation of the Sleuth 
r-band LC analysis. 



4. Spectroscopic Observations 



T-Lyrl-17236 was obs erved spectroscopica l ly wit h two instruments: The Near-Infrared 
Spectrometer (NIRSPEC; iMcLean et aD Il998l . boooh at the W. M Keck Observatory in 
Hawaii, and the Tillinghast Reflector Echelle Spectrograph (TRES; ISzentgyorgyi fc Fiiresz 
20071 ). installed on the 1.5-meter Tillinghast telescope at the F. L. Whipple Observatory in 
Arizona. 

NIRSPEC was operated using a 3-pixel slit (0.432") and an N7 blocking filter, thus 
producing a spectral resolving power of R = A/AA ~ 25,000. The duration of the exposures, 
which ranged from 420 to 900 seconds, was adjusted according to observing conditions. The 
spectra were gathered in two consecutive nods, producing a total of five NIRSPEC nod pairs. 
The nods of each pair were then subtracted one from the other, removing much of the sky 
emission. We ext r acted the spectra of both nods using the optimal extraction procedure 
outlined in iHornd (Il986l ). and then co-added the two resulting one-dimensional spectra. We 
calibrated the wavelengths of the resulting spectrum using its atmospheric telluric features, 
and then corrected for both the telluric absorption and the blaze of the spectrograph by 
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dividing this spectrum by the spectrum of an AO V- type star (HR 5511). Finally, we cross- 
correlated each spectrum with the spectrum of an M0.5V template star (GJ 182). To this 
end, we used a single NIRSPEC order (2290-2320 nm), which is within the f^-band, and 
has a scale of 0.0336 nmpixel"^ at its center. This order covers the CO 2-0 bandhead, which 
includes a rich forest of R-branch transition lines, as well as many telluric absorption features 
due to methane in the Earth's atmosphere. The advant ages offered by th is spectral region 



and the details of the instrument setup are described in lBlake et al.l (120081 ). 



TRES is a high-resolution fiber-fed optical echelle spectrograph designed to cover a 
large range of wavelengths (390-934 nm) in 51 orders. We employed the medium-size fiber 
(2.3") so as to cover the full stellar PSF, while providing a spectral resolving power of R 
^ 47,000. Following each of our three 900-1000 second exposures, the TRES data were 
read from a 4638 x 1090-pixel CCD, which we set to a 2x2 binning mode for a more rapid 
read-out. We then used a dedicated IRAF toolset to process and extract 51 spectral or- 
ders simultaneously, ultimately producing 2319 data points along each order. The IRAF 
processing of the TRES data involved merging the mosaic FITS files, removing cosmic ray 
hits, flattening fringing effects, and then extracting the orders. We wavelength-calibrated 
the TRES spectra using Thorium- Argon (ThAr) exposures, and then corrected the telluric 
absorption and spectroscopic blazing by dividing each spectrum by a TRES spectrum of a 
rapidly-rotating BOIV-type star (HR 264). Though TRES produces 51 spectral orders, we 
used only four of them, covering wavelengths of 665-720 nm (similar to the i?-band), and 
having a post-binning scale of ~0.0065 nmpixel"^. These orders contain a diverse array of 
absorption features, including those of TiO, Fe I, Ca I, Ni I, and Cr I. We limited ourselves 
to these orders because at shorter wavelengths there was insufficient flux from our red tar- 
get, while at longer wavelengths the spectra were dominated by telluric absorption features, 
produced largely by terrestrial O2 and H2O. We cross-correlated these four orders with the 
corresponding orders of an M1.5V template star (GJ 15A, also known as GX And A) and av 



eraged their cross-correlation functions. We repeated this final calculation using the IZucker 



( 120031 ) maximum-likelihood method, which reproduced our results to within a fraction of 



their uncertainties, although with slightly larger errors^ 

In total, we produced five RV measurements of each component with NIRSPEC and 
three with TRES. In all cases we were able to measure the RVs of both binary components 
by employing a cross-correlation method that transforms the spectra to Fourier-space using 



'^The lZuckeil(l200i method is more accurate than simple cross-correlation averaging for large N. However, 
because it takes the absolute value of the correlation, it loses some information and effectively increases the 
noise baseline. This increased noise will negate its advantage when combining a small number of correlations, 
as is the case in our TRES analysis {N — 4). 
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the Lomb-Scargle algorithm (jPress et al.lll992l ). This method allowed us to cross-correlate 
spectra with arbitrary sampling, without having to interpolate or resample them onto an 
equidistant grid. We then multiplied the Fourier-transformed target and template spectra, 
inverse- Fourier-transformed the product, and normalized it. Since the resulting two peaks 
in the cross-correlation functions were always well separated, we were able to fit each with 
a parabola, and thus measure their offsets and widths. The uncertainties of these RVs 
are somewhat difficult to determine with our procedures, but tests indicate that they are 
approximately 1.0 kms~^ and 1.4 kms~^ for the primary and secondary in our NIRSPEC 
spectra, and about 0.5 kms~^ and 1.2 kms"^ in our TRES spectra. These internal errors are 
adopted below in the spectroscopic analysis, but have relatively little effect on the results. 
Finally, the RVs were transformed to the barycentric frame, and the TRES RV measurements 
were further offset by —2.82 kms~^ in order to place them on the same reference frame as 
the NIRSPEC measurements, which were obtained with a different template (GJ 182). This 
offset was determined by including it as an additional free parameter in the Keplerian RV 
model (see §[5]). Once the offset was determined, we held its value fixed in all subsequent 
analyses. The final velocities are listed in Ta ble [3] and include th is offset. Note that these 
listed RVs are all relative to GJ 182, for which iMontes et al.l (120011 ) have measured the value 
+32.4 ± 1.0 kms-^ 



5. Orbital Analysis 



We began our analysis by determining the orbital period (P) and the epoch of primary 
eclipse (to), and constraining the eccentricity (e) of T-Lyrl-17236 through eclipse timing. 
The times of eclipse determined from our photometric observations listed in Table HI Since 
our data span 3.5 years, we were able to determine the period to an accuracy of 3 seconds 
(see Table [5]). To estimate the binary's eccentricity, we first measured the observed minus 
calculated [O—C) timing difference between the primary and secondary eclipses in all avail- 
able LCs, which provided an upper bound of |ecosa;| < 0.0008, where a; is the argument 
of periastron (see Figure [HD. Though u and e cannot be determined separately in this way, 
this result indicates that the orbit of T-Lyrl-17236 is likely to be circular or very nearly 
so. This conclusion is further supported by a weaker upper limit of |esinc<j| < 0.06, ob 



tainec 



1978 



thro ugh preliminary LC model fitting (see below). Theoretical estimates (jZahnlll977 
1993) of th is binary suggest a circularization timescale of tcirc — 390 Gyr (see also 



Devor et al.ll2008h . Being many times the age of the binary, this long timescale suggests 
that T-Lyrl-17236 formed in a circular orbit. However, this timescale value is an instanta- 
neous est imate for the current epoch, and is l ikely to have been significantly different in the 
past (see lZahn fc Bouchetlll989l : iMazehl l2008l . and references therein). Therefore, it is quite 
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possible that the binary circularized while it was in the pre-main sequence, however, to the 
extent that this theory is correct, it is unlikely to have circularized once settling on the main 
sequence. 

A Keplerian model was fitted to the radial velocities to determine the elements of the 
spectroscopic orbit of T-Lyrl- 17236. We assumed the eccentricity to be zero based on the 
evidence above and the lack of any indications to the contrary from preliminary spectroscopic 
solutions. The period and to were held fixed at the values determined above. We solved 
simultaneously for the velocity semi-amplitudes of the components (Ka^b) and the RV of 
their center of mass {Vy). The results are shown graphically in Figure [71 and the elements 
are listed in Table [61 The minimum masses Ma_,b sin^ i are formally determined to better 
than 2%. However, because of the small number of observations (A^ = 8), the possibility of 
systematic errors cannot be ruled out and further observations are encouraged to confirm 
the accuracy of these results. 

We then proceeded to find the remaining photometric para meters of T-Lyrl-17236. To 



this end, we analyzed the Sleuth r-band LC using JKT-EBOP (ISouthworth et al. 



2004a 



a LC mode l ing program based on the EPOB light curve generator ( Nelson fc DavisI 



1972 



Etzell Il98ll : [Popper &: Etzell Il98ll ). We assumed a circular orbit, as before, a mass ratio 
of g = 0.7692 from the spectroscopic model, and the period determined above. We solved 
simultaneously for the orbital inclination {i), the fractional radii {rA,B), the central surface 
brightness ratio of the secondary in units of the primary (J), the time of primary eclipse 
(to), and the out-of-eclipse magnitude (zero point). We estimated the uncertainties of the 
f itted parameters by eva luating the distribution generated by 1000 Monte Carlo simulations 



(ISouthworth et al.[ [20051 ). 



Because of the large photometric aperture of Sleuth, the presence of significant con- 
tamination from the light of additional stars is a distinct possibility. Unfortunately, due 
to its degeneracy with the orbital inclination and the fractional radii, we were not able to 
simultaneously determine the fractional third light of the system (Z3). We therefore sequen- 
tially refit the LC model parameters with fixed fractional third-light values ranging from 
to 0.2 (see Figure [S]). We repeated this routine with the /-band IAC80/HATNet LC as 
well, although these results were not used because of their larger uncertainties. We obtained 
an external estim ate of the third-lig ht fraction affecting the Sleuth observations using the 



an external estim ate 01 tne tnira-ug nt traction artectmg tne bieutn oDservations using tne 
USNO-B catalog jMonet et al.lboosh . which lists two dim objects within 30" of T-Lyrl-17236 



(USNO-Bl.O 1366-0314297 and 1366-0314302). Assuming that these objects are completely 
blended into T-Lyrl-17236, we expect an i?-band third-light fraction of ^3 = 0.085 ± 0.018, 
and we adopted this value for the r-band LC. Fortunately, the fitted parameters are quite 
insensitive to third light, so that the uncertainty in I3 only moderately increases their un- 



- 9 - 



certainties. No objects were listed within the smaller photometric apertures of either IAC80 
or HATNet, so we conclude that the J-band LC should have little or no third-light contam- 
ination. It is important to note that these third-light estimates assume that there are no 
further unresolved luminous objects that are blended with T-Lyrl-17236 (e.g., a hierarchical 
tertiary component). However, the divergence of the r-band and J-band solutions at higher 
third-light fractions (see Figure [8]) , and the deep primary eclipse in both the r- and /-bands 
(0.649 mag and 0.604 mag, respectively), suggest that if such unresolved objects exist, they 
are unlikely to account for more than ~0.1 of the total flux, and therefore would not bias 
the fitted results beyond the current estimated uncertainties. The final results of our LC fits 
are given in Table [51 



6. Physical Parameters 



The fundamental parameters of T-Lyrl-17236, such as their absolute masses and radii, 
were derived by combining the results of the spectroscopic analysis (Table [H]) with those 
from the photometric analysis (Table [5]). These and other physical properties are listed in 
Table [71 Our estimates of the primary and secondary component masses. Ma = 0.6795 ± 
0.0107 Mq and Mb = 0.5226 ± 0.0061 lead us t o infer spectral types of K5V and MOV, 
respectively, according to empirical tables (ICoxll2000l ). We are not able to make independent 
estimates of the effective temperatures of the stars from the data in hand. This could be 
done, for example, if we had individual color indices based on combined light values and 
light ratios in two different bands, but we can only derive a reliab le estimate of the li ght 
ratio in the r-band. The comparison with stellar evolution models by lBaraffe et al.l (119981 ) in 
§[H1 suggests primary and secondary component temperatures of approximately 4150 K and 
3700 K, respectively, although the accuracy of these values is difficult to assess. 

No trigonometric parallax is available for T-Lyrl-17236. A rough distance estimate to 
the system may be made using the JHKs brightness measurements in the 2MASS Catalog, 
collected in Table [H along with estimates of the absolute magnitudes. For these we must 
rely once again on models. The Galactic latitude of - j-16.8° suggests the p ossibility of some 
interstellar extinction. From the reddening maps of ISchlegel et al.l (Il998l ) we infer E{B — 
V) ~ 0.07 in the direction of the object (total reddening), which corres ponds to e xtinctions of 
A{J) ~ 0.061, A{H) ~ 0.038, and A{K) ~ 0.011, assuming Ry = 3.1 jcoxlboooh . Under the 
further assumption that this extinction applies to T-Lyrl-17236, we derive a mean distance 
of 230 ifc 20 p c, after conversion of the near-infrared magnitude s in the CIT syste m from the 
Baraffe et all Jl998h models to the 2MASS system, following [CarEenterl koOlh . With the 
proper motion components from the USNO-B Catalog listed in Table [2l the center-of-mass 
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veloc ity V^y from the spectroscopic solution corrected for the velocity of GJ 182 (IMontes et al. 



20011 ). and the distance above, we infer space velocity components in the Galactic frame of 



{U,V,W) ~ (+41, +21, +2) kms ^, where U points in the direction of the Galactic center. 

Because of the relevance of the rotational velocities of the stars for the interpretation 
of the chromospheric activity results of §[7l we have made an effort here to measure the 
rotational broadening of both components from the widths of the cross-correlation functions 
derived from our TRES spectra. We rely on the fact that to first order, the width of a cross- 
correlation peak is approximately equal to the quadrature sum of the line broadening of the 
two spectra. We began our estimation procedure by finding the effective resolution of the 
instrument (cTj) in the four TRES orders we used. This was done by auto-correlating a TRES 
ThAr spectrum that was taken just before the second T-Lyr 1-17236 observation. We found 
that the four orders produced peaks with an average FWHM of 8.90 ± 0.17 kms~^. Thus, 
assuming that the intrinsic widths of the ThAr emission lines are negligible compared to 
the instrumental resolution, we found that cxj = 6.29 ± 0.12 kms~^. This value corresponds 
to a spectral resolving power of R = 47,630 ± 930, which is consistent with the TRES 
specifications. Next, we determined the intrinsic spectral line broadening of the template 
star, GJ 15A (cij). We auto-correlated the template spectrum and found that it produced 
peaks with an average FWHM of 9.7 ± 1.4 kms'^ This value should be equal to V2{af + 
^2^1/2^ from which we infer th at at = 2.7 ± 2 .5 kms ~^. Note that this result is well within 



the upper bound provided by iDelfosse et al.l ( 119981 ). following their non-detection of any 



rotational broadening in GJ 15A. Using this information, we can now find the intrinsic 
spectral line broadening of the T-Lyr 1-17236 components {o'a,b)- The average FWHM of the 
primary and secondary peaks, resulting from the cross- correlation of each observed spectrum 
of T-Lyrl-17236 against the template, were measured to be 12.6 ± 2.0 kms~^ and 12.0 ± 
2.4 kms""*^, respectively. These widths are expected to be equal to + o"^) + (crf + cr^^)]^/^, 
from which we calculate that a a = 8.4 ± 3.0 kms~^ and as = 7.6 ± 3.8 kms^^. 

The rotational profile FWHM expected for a homogeneous stellar disk is y/Svsinir, 
where v is the star's equatorial rotational velocity, and ir is the inclination of its rotational 
axis. Stellar limb dar kening, how ever, will narrow the rotational profile, thus decreasing 



the observed F WHM (|Gravl 119921) . Adopting the i?-band PHOENIX linear limb darkening 
coefficients from Claret ni998h . we find that the expected FWHM values for the primary 



and secondary components of T-Lyrl-17236 are, respectively, 1.495 f sin and 1.499 f sin v. 
Using these results we can set upper bounds to the components' v sin These upper bounds 
represent the limiting case whereby the spectral line broadening is due entirely to stellar 
rotation, and we neglect all other line broadening mechanisms, such as microturbulence and 
the Zeeman effect. We thus determine the maximum rotational velocities of the T-Lyrl-17236 
primary and secondary components to be vsinir = 5.6 ± 2.0 kms~^ and 5.1 ± 2.3 kms~^. 
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respectively. 



An estimate of the timescale for tidal synchronizat ion of the st ars' rotation with their 
orbital motion may be obtained from theor y following IZahnI (119771 ). and assuming simple 
power-law mass-radius-luminosity relations (jCoxl 120001 ). Thus, for stars less massive than 
1.3 Mq, 



''sync — 



0.00672 Myr (A;2/0.005) 



-\-\l + qf{P/d^yf{M/M^) 



-4.82 



1^ 



where k2 is determined by the structure a nd dynamic s of the star and can be obtained by 
interpolating published theoretical tables (jZahnlll994l ). This calculation leads to timescales 
of tsync — 0.56 Gyr and 1.02 Gyr for the primary and secondary components of T-Lyrl-17236, 
respectively, which are much shorter than the circularization timescale determined in §[5l We 
note that similar to the circularization timescale, the synchronization timescales estimated 
above are the current instantaneous values, and are likely to have changed over time. The 
age of the system is undetermined (see §[8]), but assuming its age is at least a few Gyr, as is 
typical for field stars, it would not be surprising if tidal forces between the components had 
already synchronized their rotations. This is illustrated in Figure [H where T-Lyrl-17236 
is shown along with the other systems in Table [T] and with curves representing theoretical 
estimates of the synchronization timescale as a function of orbital period. 

If we assume that the components are indeed rotationally synchronized, we can compute 
their rotational velocities more accurately using va,b = '^t^Ra,b / P- We thus derive synchro- 
nized velocities of (f sinv)sync = 3.81 ± 0.26 kms~^ and 3.15 ± 0.31 kms~^ for the primary 
and secondary components, respectively. These values are slightly below but still consistent 
with the maximum rotational velocities measured above. Thus, observational evidence sug- 
gests that the stars' rotations may well be synchronized with their orbital motion, although 
more precise measurements would be needed to confirm this. Our conclusion from this cal- 
culation is that regardless of whether we assume that the components of T-Lyrl-17236 are 
synchronized, their rotational velocities do not appear to be large. 



7. Chromospheric Activity 

Our absolute mass and radius determinations for T-Lyrl-17236 offer the possibility of 
testing stellar evolution models in the lower main sequence, and in particular testing the idea 
that the discrepancies noted in §[T] are related to chromospheric activity and the associated 
magnetic fields in systems where the components are rotating relatively rapidly. Thus, es- 
tablishing the level of the activity in the system presented here is of considerable importance. 
We have shown in §|6]that the relatively long period of T-Lyrl-17236 (P ~ 8.429441 days) 
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implies that even if the components are synchronized, their rotational velocities are slow, 
and therefore are not expected to induce a great deal of chromospheric activity. However, 
demonstrating that the stars are indeed inactive requires more direct evidence, given that 
some stars of similar masses as these are st i ll fou nd to be quite active at rotation periods 



as long as 8 days (see, e.g., iPizzolato et al.l 120031 ). We present here the constraints avail 



able on the surface activity of T-Lyrl-17236 from its X-ray emission, optical variations, and 
spectroscopic indicators. 



The present system has no entry in the ROSAT Faint Source Catalog ( IVoges et al.lll999l ). 
suggesting the X-ray luminosity, usually associated with activity, is not strong. Examination 
of the original ROSAT archive images leads to a conservative upper limit to the X-ray flux 
of 6.71 X 10~^^ ergcm"^ s~^ in the energy range 0.1-2.4 keV, and together with information 
from Table [TJ we infer an upper limit for the ratio of the X-ray to bolometric luminosity of 
\ogLxl L\,o\ < —3.13 . Values for the four best studied cases of CM Dra, YY Gem, CU Cnc, 
and GU Boo, which are all very active, are respectively —3.15, —2.88, —3.02, and —2.90 



see 



Lopez- Morales! 120071 ) . These are at the level of our limit or higher, although we do not 



consider this evidence conclusive. 

There are no detectable variations in the r-band light curve out of eclipse, within the 
uncertainties. Such variations would be expected from activity-related surface features show- 
ing significant contrast with the photospheres. We estimate an upper limit of ~0.01 mag 
in r for the night-to-night variations (see Figure 3). Because the secondary components is 
significantly dimmer, it has a weaker variability upper limit of ~0.09 mag. We note, how- 
ever, that this evidence for inactivity is not conclusive either, since the observed photometric 
variations can depend significantly on the distribution of spots on the surface. 

A number of spectroscopic activity indicators (the Ca II H and K lines. Ha, etc.) 
should in principle allow a more direct assessment of the activity level in T-Lyrl-17236. 
Unfortunately, however, the quality of our spectroscopic material in the optical makes this 
difficult. The flux in the blue for this very red system is too low to distinguish the Ca II H 
and K lines, and even at Ha the noise is considerable (typical signal-to-noise ratios at this 
wavelength are ~12 pixel"^). Two of the three TRES spectra show the Ha line in absorption, 
and the other appears to show Ha in emission. This suggests some degree of chromospheric 
activity, although perhaps not at such a high level as to sustain the emission at all times, as 
is seen in other stars. H/3 appears to be in absorption in all three TRES spectra. 

Clearly more spectra with higher signal-to-noise ratios are needed to better characterize 
the level of activity, but from the sum of the evidence above it would not appear that the 
activity in T-Lyrl-17236 is as high as in other low-mass eclipsing binaries studied previously, 
thus more closely aligning it with the assumptions of current standard stellar models. The 
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system may therefore constitute a useful test case for confirming or refuting the magnetic 
disruption hypothesis (see §[1]), which predicts that the absolute properties of its slowly 
rotating components should match the theoretical models of convective stars. 



8. Comparison with Models and Conclusions 



A comparison with solar-metallicity models by iBaraffe et al.l (Il998l ) for a mixing length 
parameter of oml = 1-0 is presented in Figure [91 Our mass and radius determinations for 
T-Lyrl-17236 (see Table [7]) are shown along with those of the low-mass systems listed in Ta- 
ble m The location of the models in this diagram depends only slightly on age because these 
stars evolve very slowly. The age of T-Lyrl-17236 is difficult to establish independently. The 
space motions derived in §|6]do not associate the system with any known moving group, and 
are quite typical of the thin disk. Thus, all we can say is that it is not likely to be very 
old. We display in Figure M two models for ages of 1 Gyr and 10 Gyr, which likely bracket 
the true age of T-Lyrl-17236. Within the errors, our measurements for the two components 
are consistent with the models, which would in principle support the magnetic disruption 
hypothesis. Unfortunately, however, the uncertainties in the radius measurements (~7% 
and ~10%) are still large enough that our statement cannot be made more conclusive. Fur- 
ther follow-up observations, especially rapid-cadence and precise photometric measurements 
during multiple eclipses, should significantly reduce the uncertainties in the radii and thus 
provide far stronger constraints on the theoretical models of low-mass stars. Additionally, 
higher-quality spectroscopic observations than ours are needed to confirm that the level of 
chromospheric activity in the system is relatively low. If after such observations, the masses 
and radii of the T-Lyrl-17236 components remain consistent with the stellar models, then 
the magnetic disruption hypothesis will be strengthened. However, if further observations 
find that the components of T-Lyrl-17236 are larger than predicted by current stellar mod- 
els, as is the case with most other similar systems investigated in sufficient detail, then this 
will provide evidence that additional mechanisms n eed to be included in the models of the 



structure of low-mass main sequence stars (see, e.g.. IChabrier et al.ll2007l ). 



It is important to note he re that T-Lyrl-1723 6 falls within the field of view of the up 



coming NASA Kepler Mission (IBorucki et al.ll2003[ ). The Kepler Mission will not return data 



for all stars within its field of view; rather, the targets will be selected by the Kepler team. 
We see at least two reasons why such monitoring of T-Lyrl-17236 would be of significant 
value. First, the data would greatly refine the estimates of the physical parameters of the 
component stars and may permit a search for their asteroseismological modes. Second, the 
data would enable a search for transits of exoplanets, which are expected to orbit in the 
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same plane as that defined by the stellar orbits. 

Finally, we note that our findings in this paper confirm the accuracy of the MECI 
algorithm (see Figure [2]), which can be further used to find additional long-period low-mass 
EBs, and indeed a variety of other interesting targets. We have shown in a recent paper 



(jPevor et al.ll2008l ) how this can be done with comparable ease by systematically searching 
the ever-growing body of LC survey datasets. We hope that this new approach for locating 
rare EBs will motivate additional studies of these vast, largely untapped datasets, which 
likely harbor a wealth of information on the formation, structure, dynamics, and evolution 
of stars. 
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Table 1. Periods of well characterized main sequence EBs with both component masses 





below O.S 




Name 


Period [days] 


Citation 


OGLE BW5 VaS" 


0.198 


Maceroni & Montalban (2004) 


RR Caeli'' 


0.304 


Maxtcd et al. f2007') 


NSVS01031772 


0.368 


Lopez-Morales et al. ("2006) 


SDSS-MEB-1 


0.407 


Blake et al. (2007) 


GU Boo 


0.489 


^opez-Moralcs &^^^as ^00^ 


2MASS J04463285+1901432 


0.619 


^ebb ot al. (2006^ 


YY Gem 


0.814 


^ron (1952); Torres & Ribas (2002) 


T-HerO-07621 


1.121 


Creevev et al. (2005) 


CM Dra 


1.268 


Lacv (1977b): Metcalfe et al. (1996) 


UNSW-TR-2 


2.117 


Youne et al. (2006) 


2MASS J01542930+0053266 


2.639 


Becker ot al. (2008) 


CU One 


2.771 


Delfosse et al. (1999): Ribas (2003) 



'^This binary might not be detached, as its components seem to be undergoing significant 
mutual heating and tidal interactions due to their proximity (a — 1.355 i 0.066-Rq). 

^This is an unusual case of an EB containing a white-dwarf (primary) and an M-dwarf 
(secondary). As such, the primary component is likely to have transferred mass to the 
secondary component, and perhaps even enveloped it during the red-giant phase of its 
evolution. 



Table 2. Catalog information for T-Lyrl-17236 



Source Catalog 


Parameter 


Value 


2MASS'' 


Q {J2000) 


19:07:16.621 


2MASS 


S (J2000) 


-1-46:39:53.21 


usno-b'' 


B mag 


16.11 ± 0.2 


GSC2.3" 


V mag 


14.37 ± 0.28 


USNO-B 


R mag 


14.41 ± 0.2 


CMC14'^ 


r' mag 


14.073 ± 0.029 


2MASS 


J mag 


12.019 ± 0.015 


2MASS 


H mag 


11.399 ± 0.015 


2MASS 


Ks mag 


11.235 ± 0.015 


USNO-B 


(masyr~"^) 


-2 ± 3 


USNO-B 


MS (masyr^^) 


-28 ± 2 


2MASS 


identification 


19071662-1-4639532 


CMC14 


identification 


190716.6+463953 


GSC2.3 


identification 


N2EH033540 


USNO-B 


identification 


1366-0314305 



^Two Micron All Sky Survey catalog jSkrutskie et al.| 
120061 '). 

U.S. Naval Observatory photographic sky survey 
jMonet et al.ll2003l) . 

'^Guide Star Catalog, version 2.3.2 jMorrison et al.ir200l|') . 

'^Carlsberg Meridian Catalog 14 jEvans et al.ll2002ll . 
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Table 3. Radial velocity measurements for T-Lyrl-17236 in the barycentric frame, relative 

to GJ 182 



1 /^Tx^\ Prinrary RV Secondary RV Exposure Time rr^ . , t , 

Epoch (BJD) , _i , , _i , . ■. Template Instrument 

(km s ^ ) (km s ^ ) (sec) 



2453927.9400 


-2.87 


-45.24 


480 


GJ 182 


NIRSPEC 


2453930.9258 


-68.09 


38.86 


900 


GJ 182 


NIRSPEC 


2453946.8846 


-64.26 


36.63 


600 


GJ 182 


NIRSPEC 


2453948.9100 


-43.46 


7.03 


420 


GJ 182 


NIRSPEC 


2464312.7985 


7.66 


-67.68 


480 


GJ 182 


NIRSPEC 


2454372.6179 


23.99 


-80.14 


900 


GJ 16A 


TRES 


2454377.6382 


-68.03 


40.10 


1000 


GJ 16A 


TRES 


2454377.6624 


-67.97 


39.73 


1000 


GJ 16A 


TRES 



Table 4. Eclipse timings measured for T-Lyrl-17236 



Eclipse Type Epoch (HJD) O-C [sec] Data Source 



Primary 


2463162, 


.96121 


900+232 


HATNet 


Secondary 


2463167, 


.17693 


—646+^*8* 


HATNet 


Primary 


2463169, 


.82009 


48+126 
*°-131 


HATNet 


Primary 


2463186, 


.67897 


„,7+214 
^=•^'-221 


HATNet 


Secondary 


2463190. 


.89369 


,,1+431 
— 423 


HATNet 


Primary 


2463196, 


.10841 




HATNet 


Secondary 


2463207, 


.76268 


„„,,+642 
^=^0-648 


HATNet 


Secondary 


2453544, 


,93022 


^'^^-332 


Sleuth 


Secondary 


2453561 


,78910 


312 + °^ 


Sleuth + PSST 


.Secondary 


2453578 


,64798 


-,, + 206 

''-^"•-208 


Sleuth 


Primary 


2453582 


,86270 




Sleuth 


Primary 


2453599, 


,72158 


941^4 


Sleuth 


Secondary 


2453603, 


,93630 


1047t-t 


Sleuth + PSST 


Primary 


2453616, 


,58046 


■"-175 


Sleuth 


Primary 


245;-i8fil 


I):i425 


23_^+28n 


PSST 


Prinialx 


2 ■i."'. -1-1 17 






lACsd 



Table 5. Photometric parameters of T-Lyrl-17236 



Parameter 


Symbol 


Value 


Period (days) 


P 


8.429441 ± 0.000033 


Epoch of eclipse (HJD) 


to 


2453700.87725 ± 0.00041 


Primary fractional radius 


rA 


0.0342 ± 0.0023 


Secondary fractional radius 




0.0283 ± 0,0028 


Orbital inclination [deg] 


i 


89.02 ± 0.26 


Eccentricity 


e 


0.0 (fixed) 


Sum of fractional radii 


rA + rB 


0.06266 ± 0.00096 


Ratio of radii (Rb/Ra) 


k 


0.83 ± 0.16 


Light ratio (r-band) 


^bI^a 


0.173 ± 0.073 


Surface brightness ratio (r-band) 


Jb/ Ja 


0.2526 ± 0.0099 
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Table 6. Spectroscopic parameters of T-Lyrl-17236 



Parameter 


Symbol 


Value 


Primary radial velocity semi-amplitude (kms "'') 


Ka 


48.36 ± 0.23 


Secondary radial velocity scmi-amplitudc {kms""'") 


Kb 


62.86 ± 0.46 


Baryccntric radial velocity, relative to GJ 182^ {kms "'" ) 


v-i 


-21.01 ± 0.18 


Binary separation with projection factor (-Rq) 


a sin i 


18.526 ± 0.083 


Primary mass with projection factor (Mq) 


Ma sin^ i 


0.6792 ± 0.0107 


Secondary mass with projection factor (Mq) 


Mb sin"^ i 


0.5224 ± 0.0061 


Mass ratio {Af^/M^) 


1 


0.7692 ± 0.0069 



iMontea et al.l )200lD list the radial velocity of GJ 182 as +32.4 ± 1.0 kms^^. 



Table 7. System parameters of T-Lyrl-17236 



Parameter Symbol Component A Component B 



Mass (Mq) 


M 


0.6795 ± 0.0107 


0.5226 ± 0.0061 


Radius (Hq) 


R 


0.634 ± 0.043 


0.525 ± 0.052 


Log surface gravity (cgs) 


log a 


4.666 ± 0.059 


4.718 ± 0.086 


Semimajor axis (10^ km) 


a 


5.606 ± 0.027 


7.288 ± 0.053 


Maximum rotational velocity^ (kms~"^) 


V sin V 


5.6 ± 2.0 


5.1 ± 2.3 


Synchronized rotational velocity^ (kms~^) 


(^v sin ii')sync 


3.81 ± 0.26 


3.15 ± 0.31 


Absolute visual magnitude^ (mag) 


Mv 


8.03 


9.67 


Bolometric luminosity^ (^0) 


L 


0.110 


0.039 


Effective temperature^ (K) 




4150 


3700 


Distance^ (pc) 


D 


230 ± 20 



^Sce description in ^6] 

Inferred using stellar evolution mo dels by iBaraffe et al.l )1998D assuming solar metallicity and an age 
of 2.5 Gyr. 
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2 4 6 8 

Orbital Period [days] 



Fig. 1. — The predicted synchronization timescales due to turbulent dissipation (jZahnlll977l . 



19941 ) for well characterized low-mass EBs from Table [H The lines trace constant synchro- 
nization timescales of binary components for which q = 1 (see ^ for further details on this 
calculation). The black circles indicate primary components and the grey circles indicate 
secondary components. Note that in some cases the primary and secondary symbols nearly 
overlap. 
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— ^ ^ — 








0.4 0.6 

Primary Component Mass [M©] 



Fig. 2. — The mas s- mass likelihood plot for T-Lyrl- 17236 created with MECI, using the 



Baraffe et al.l (119981 ) isochrones for an age of 2.5 Gyr. This analysis incorporated the r-band 
LC and the 2MASS colors of the target. The contour lines indicat e the weighted reduced 
chi-squ ared values of each component mass pairing, using w = 10 (jPevor fc Charbonneau 
2006bl ). The white point indicates the our final mass estimate from this paper, and the white 
square approximates our current mass uncertainties. 
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Fig. 3. — Phased light curves of T-Lyrl-17236. The top curve is from the IAC80 (black) and 
HATNet (grey) telescopes, both of which observed in J-band. Note the tight cluster of IAC80 
observations near phase 0.7; these points determine the IAC80 LC zero point. The bottom 
curve is from the Sleuth (black) and PSST (grey) telescopes, which observe, respectively, in 
the r-band and i?-band. The secondary eclipse is about twice as deep in the J-band as it 
is in the r- or i?-bands, indicating that the secondary component is significantly redder and 
therefore cooler than the primary. 
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Fig. 4. — Enlargement of the eclipse phases in the LC of T-Lyrl-17236, as recorded by the 
Sleuth (black) and PSST (grey) telescopes (r-band and i?-band, respectively). The solid line 
shows the best-fit JKT-EBOP model, for which the residuals are displayed at the bottom. 
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Fig. 5. — Enlargement of the eclipse phases in the LC of T-Lyrl-17236, as recorded by the 
IAC80 (black) and the HATNet (grey) telescopes (/-band). The solid line shows the best-fit 
JKT-EBOP model, for which the residuals are displayed at the bottom. 



-26- 




HJD - 2453000 (P = 8.429 days) 



Fig. 6. — Eclipse timing (0—C) measurements of T-Lyrl-17236. The solid triangles indicate 
primary eclipses, and the starred symbols indicate secondary eclipses. The large error bars 
are generally due to eclipses that are constrained by only a few observations, or for which only 
the ingress or egress was observed. The cluster of points at the very left (HJD < 2,453,300) 
are measurements from HATNet, the single data point at HJD 2,454,417 is from the IAC80, 
and the remaining data are from Sleuth and PSST. The two parallel dashed hnes indicate 
the expected O — C location of the primary (bottom) and secondary (top) eclipses, in the 
best- fit eccentric model (|ecosa;| ~ 0.0005). This eccentric model provides only a very 
small improvement in the fit compared to the circular model (F-test: ')d,^circl 'xl,ecc — 1-29, 
indicating a p ~ 0.33 significance). 
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Fig. 7. — Radial velocity measurements of T-Lyrl-17236, relative to GJ 182, shown as a 
function of orbital phase. The velocities of the primary component are represented with 
squares, and those of the secondary with circles. The filled symbols correspond to data 
taken with NIRSPEC, and the open symbols represent TRES measurements. Residuals 
from the model fit are shown below for the primary and secondary components. 
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Fig. 8. — The JKT-EBOP parameter fits over a range of values for the third hght fraction. 
The panels show the best-fit values and uncertainties for: (A) The fractional radii of the 
primary (r^), and (B) secondary (r^) components; (C) The sum of the fractional radii 
{ta + i^b)) and (D) the radius ratio {k = tb/ta)-^ (E) The binary orbital inclination (i), and 
(F) the central surface brightness ratio (J, secondary over primary). Note that in contrast 
to other panels, panel (F) shows distinct values for the /- and r-band LCs. This is expected 
since the two components have different colors, and therefore different relative fluxes through 
different fllters. In all cases the estimated third light fractions for the r-band and the J-band 
LCs are indicated by boxes. 
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Fig. 9. — Mass-radius diagram for T-Lyr 1-17236 and other low-mass ecl ipsing binaries under 
0.8 Mq from Tabled! Theoretical isochrones for solar metallicity from iBaraffe et al.l (119981 ) 
are shown for ages of 1 and 10 Gyr. The components of T-Lyrl-17236 are indicated with 
arrows. Most of these binary components (particularly those with smaller uncertainties) 
display a systematic offset in which their measured radii are larger than predicted from 
models. 



